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1. Introduction and the model 

Many optical media are based on Bragg gratings (BGs), which are built as modulations of the local refractive index 
with the spatial period close to the half-wavelength of the electromagnetic field. Fiber BGs, which operate in the 
temporal domain, are used in numerous applications, and are also a subject of fundamental studies M- The 
interplay of the BG with the nonlinearity of the fiber material gives rise to gap solitons S-Si, which were originally 
predicted as analytical solutions and then created in experiments [l(3]-[lfj. Spatial BG solitons were predicted 

too, using periodic gratings written on a planar waveguide [lTf-fl^] or in a photonic crystal [20]. Solitons can also be 
supported by a combination of the BG with the second-harmonic-generating 21|-[24| or resonant (two-level) (25l. (26| 
nonlinearity. 

The variety of modes supported by the BGs and their potential applications can be expanded by using various 
superstructures created on top of regular gratings. The technology of the creation of such complex gratings in optical 
fibers became available long ago [27]-In particular, supergratin gs g ive rise to additional spectral bandgaps, in 
addition to the central one induced by the underlying uniform BG [30|_. Theoretical models of the supergratings, 
incorporating the nonlinearity of the material, were developed in detail [30|-j32]. 

Another example of the superstructure is known in the form of the Moire pattern, with a sinusoidal modulation 
imposed on the periodic variation of the refractive index building the ordinary BG. It gives rise to a narrow transmission 
band in the middle_of the central gap (resembling the effect of the electromagnetically-induced transparency), which 
was proposed (H, [34] and implemented experimentally (35| for the retardation of light in the BG. 

It is also possible to design BG superstructures patterns in a “semi-discrete” form, i.e., as a uniform nonlinear 
waveguide with periodically inserted narrow Bragg reflectors [36j , and as a more general scheme, which makes use of 
a periodic modulation of the BG chirp and local Bragg reflectivity [37j . On the other hand, complex grating patterns 
can be readily written in photonic crystals and photonic-crystal fibers fdBMdi)} . 

Another relevant setting is a uniform nonlinear BG with one or several local defects (rather than defect lattices) 
4i|-[45j. In particular, solitons pinned to a defect represented by a ^-functional perturbation of the local refractive 
index or Bragg reflectivity can be found in an exact analytical form [42] . The delta-function approximates a BG 
with a narrow embedded layer carrying large reflectivity of the opposite sign. Interesting results are also produced 
by models which include symmetric pairs of defects: repulsive ones form cavities trapping solitons which perform 
shuttle oscillations, while attractive pairs give rise to pinned states featuring spontaneous symmetry breaking [47j . 
The interaction of moving BG solitons with weak local defects can be described by means of the perturbation theory 
which was developed in detail in Refs. [48]-[50|. 

The nonlinear light propagation in the “apodized” BG, i.e., one with a spatially inhomogeneous Bragg reflectivity, 
obeys the system of coupled-mode equations for amplitudes u and v of the right- and left-traveling waves [3l|, l48B50j : 


i 


du 

dt 


.du 


k{x)v 



u = 0, 


.dv .dv 
l ~dt~ l 'd^ 


K*(x)u 



v = 0, 


( 1 ) 


where t and x are the time and coordinate normalized so that the group velocities of the waves are ±1, ratio 2 : 1 of 
the cross- and self-phase-modulation coefficients in the cubic terms is a standard feature of the Kerr effect [HJ, and 
k(x) is the inhomogeneous Bragg reflectivity, which may be complex, taking into regard a phase shift between the 
incident and reflected waves (the asterisk stands for the complex conjugate). Replacing t by propagation distance z, 
the same coupled-mode system may be realized in the spatial domain, i.e., for planar waveguides instead of fibers. 

A local finite-width defect corresponding to 


re(x) = 1 — 6 sech (kx ), (2) 

which describes a local depression of the Bragg reflectivity, was studied in Ref. [44[. Because such a defect attracts 
the wave fields, it can support localized pinned modes, including defect modes (DMs) in the linearized version of the 
model. The DMs and their nonlinear extension were constructed in Ref. 144^ by means of numerical methods, and the 
stability of the nonlinear pinned modes was investigated. It was found that the instability is caused by resonances of 
perturbations with the continuous spectrum. 

The aim of the present work is to introduce a similar but simpler and, in a certain sense, more general, configuration, 
with a gapless layer placed between two semi-infinite gratings. We define two versions of this system. In one, the 
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semi-infinite BGs are mutually symmetric, with a possible phase shift (a) between them: 


k(x) 


1, at x > 0, 

0, at — L < x < 0, 
exp (ia) , at x < —L, 


( 3 ) 


where |re| = 1 for k ^ 0 may be fixed by obvious rescaling, provided that reflectivity |k| takes equal absolute values 
in both semi-infinite gratings. The second system is composed of asymmetric BGs, with unequal reflectivities, k = 1 
and k = cos a (here we set 0 < a < 7t/2): 


{ 1, at x > 0, 

0, at — L < x < 0, (4) 

cos a, at x < —L. 

The presence of parameter a in either version, based on Eq. © or ®> is an essential generalization in comparison 
with the model based on Eq. ©. A still more general configuration, with a phase shift between asymmetric BGs, 
is possible too, but such system with the double complexity is beyond the scope of the present paper. Another 
advantage of the defects with the simple shapes defined by Eqs. © and © is that it is much easier to realize them 
in the experiment. Further, we demonstrate below that both linear DMs and gap solitons pinned by the gapless 
layer can be found in an analytical form, which is not possible in the case of the more sophisticated profile ©. In 
particular, we conclude that the nonlinear extension of the DMs fully merges into families of the pinned solitons. 

It is necessary to mention that Eq. © neglects the ordinary dispersion or diffraction in the gapless layer, in 
the temporal and spatial implementations of the BG, respectively. This assumption is relevant if the respective 
dispersion/diffraction length, corresponding to L , is much larger than the characteristic soliton-formation length, z so i. 
Note that, in physical units, L = 1 (this scale plays a dominant role below) corresponds to length ~ 1 mm, in the 
spatial and temporal domains alike, and, in either case, z so i takes values between 1 mm and 1 cm Q-S, 00. 
With the carrier wavelength — 1.5 /im, L ~ 1 mm corresponds to diffraction length ~ 10 cm z so i. In the temporal 
domain, the disparity between the characteristic dispersion length and z so i is sill larger. 

The diffraction in the gapless region should be taken into regard in the limit of L —> oo, which actually corresponds 
to an interface (surface) separating a semi-infinite BG and the uniform medium. The interaction of gap surface 
solitons with such an interface was studied before (Td ITbl . 

The subsequent presentation is structured as follows. Analytical results are presented in Section II. They include 
exact solutions for linear DMs in both models © and © , a perturbative treatment of its weakly nonlinear extension 
(explicitly presented for the symmetric system, with a = 0), which is well corroborated by numerical results, and an 
implicit analytical solution for solitons pinned to the gapless layer. Section III reports basic numerical results for the 
nonlinear asymmetric system, based on Eq. a, including its symmetric limit, with a = 0. The numerical results 
display existence and stability areas for the pinned solitons in the underlying parameter space, and simulations of 
collisions of moving solitons with the gapless layer (for another type of the defect, collisions were simulated in Refs. 
fill and [0). In particular, the collision may generate a trapped soliton from a part of the energy of the incident one, 
which can be used for the creation of standing optical solitons (thus far, the best experimental result for slow-light 
BG solitons was the creation of one moving at velocity ~ 0.16 speed of light in vacuum 0 ). 

2. Analytical results 
2.A. Stationary equations 

Stationary solutions of Eq. © with frequency ui are looked for as 

{u (x , t) ,v(x,t)} = e~ lu,t {U(x),V (a:)} , (5) 

where complex functions U and V satisfy the system of ordinary differential equations: 


.dU 

+i— — I - coU + k(x)V + 
ax 


\v\ 2 + l -\u \ 2 


17 = 0, 


dV 


+ ouV + n*(x)U + 


(W+^i© 2 ) 


E = 0, 


which are compatible with the constraint reducing the two stationary fields to one, 


( 6 ) 


U* (x) = -V(x). 


( 7 ) 
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The stationary solutions are characterized by the total energy, alias norm (or total power, in terms of the spatial- 
domain model) and Hamiltonian, 


E = 


+OO 


\U(x)\ 2 + \V(x)\‘ 


dx, 


( 8 ) 


H = 



u d -^-v dV 


dx 


dx 


+ c.c. 

*+oo 


(|t /| 4 + |H | 4 + |C/| 2 |H| 2 ))^ + •^Bragg 


/ +oo 

[k(x)U*V + c.c.] dx, 

-oo 


(9) 

( 10 ) 


where c.c., as well as *, stand for the complex conjugate. 

The well-known analytical solution for the quiescent BG solitons in the uniform grating, with L = 0 and constant 
k = |k| e la , corresponds to general expressions © and ©, with [G, §j, @] 


tj ,_ 1 / 2 i.'| c W 2 sinf? 

so V 3 cosh (|k| a;sin 6 — iO/2) ’ 

(ID 

cu so i = |k| cos 9. 

(12) 

Here, 6 is an intrinsic parameter of the soliton family which takes values 0 < 9 < ir, 
frequency (fl2l) belonging to the bandgap of the linearized version of Eqs. ©, 

with the respective soliton’s 

— |k| < u < + \k\. 

(13) 

The total energy © of this exact solution is 


£sd = (8/3)0 

(14) 

(note that E so i does not explicitly depend on |k|). It is known that the exact solitons are stable, in the uniform 
medium, at @]-[54j, Q] 

9 < 0max W 1.01 X 7t/2, 

(15) 

i.e., in the interval of frequencies 


(w so i) cr « -0.016 < w so i < 1, 

(16) 


cf. the full bandgap (fl3l) . 

In the asymmetric system with reflectivity profile ©, the full bandgap is reduced, per Eq. m, from w 2 < 1 to 
w 2 < cos 2 a. In the interval between the bandgaps of the two semi-infinite gratings, i.e., at 

cos 2 a < w 2 < 1, (17) 

the asymmetric system gives rise to states which are localized at x > 0 but delocalized at x < —L, see below. 

Exact solution hid features an obvious symmetry, {[/ so i(— x), y so i(—x)} = {f/* ol (x), T7*i(x)}. Similarly, localized 
modes in the system based on Eqs. © and © with a = 0 are symmetric with respect to the midpoint of the gapless 
layer, x = —L/ 2: 


{U(-x), V(-x)} = {U*(x - L ), V*(x - L)} . (18) 

The presence of the phase shift (a) in Eq. © between the semi-infinite gratings affects this relation (see below), but 
the corresponding intensity profile remains symmetric: 

|U(-:r)| 2 + \V(-x)\ 2 = \U(x - L )| 2 + \V(x - L )| 2 . 


(19) 










5 


2.B. The linear defect mode 
2.B.I. The symmetric system 

The linearized version of Eq. ©, subject to reduction 0, reduces to the single linear equation, 

i ——K uiU — k(x)U* = 0. (20) 

This equation, with k(x) taken as per Eq. © for the symmetric system, including the phase shift, gives rise to an 
exact composite (three-layer) DM solution with eigenfrequency aj e iger3 


U(x) = eU\(x) = e 


exp 


exp ^ COS (cCejgen 

S D M exp (iWeigen (x + ■§)) 
(k^eigen) 


- a / ! - ^ei S en X » at * > °- 


4 

at 


— 7} cos 


1 - W eigen (* + L) 


— L < x < 0, 

, at x < —L, 


( 21 ) 


where e is an infinitesimal amplitude which underlies the linearization. Relative amplitude Rdm and the eigenfrequency 
[it always falls into bandgap (fTH ] are determined by the condition of the continuity of U{x) at x = 0 and x = —L: 


Bdm = exp (ia/ 4 + inn ), 


( 22 ) 


k^eigen-^ H - (ck/2) — COS (^eigen) “1“ 27T77-, (23) 

where n is an arbitrary integer. Note that the DM given by Eq. m is fully symmetric about the midpoint, as 
per definition (11811 . at a = 0, and satisfies symmetry relation (1191) at any a. The respective integral energy © is 

Edm = e: 2 [l + 1 /- Wg igen j ■ 

For a weak gapless layer, i.e., L —> 0 and a —> 0, Eq. (l23l) for the fundamental DM (with n = 0) yields 

w eig en~l-(l/8)(a + 2L) 2 , (24) 

with solution m degenerating into a constant at L = a = 0. On the other hand, at a —> n, Eq. (l23l) with n = 0 
yields a vanishingly small eigenfrequency, 

_ 1 7T- a , 0 s 

Weigen ~ 2 \ L' 

Further, the largest possible eigenfrequency, w e i ge n = 1, is attained by the solution to Eq. (l23l) for the fundamental 
DM (n = 0) at a (w e i g en = 1) = —2 L, and close to this point the solution takes the same form as given by Eq. (1241) . 
The fundamental DM does not exist at 

a < a m i n = -2 L. (26) 

The relation between the eigenfrequency and a, produced by Eq. (1?51) with n = 0, is displayed, for L = 0.5 
and L = 2.0, in Fig. [Tj [oj < 0 corresponds to the DM subject to constraint U(x) = +V*(x), instead of Eq. J7])], 
which demonstrates that the relation keeps its qualitative shape at all relevant values of the gapless-layer’s width. In 
particular, in the limit of L —> 0, which corresponds to the junction of two semi-infinite gratings with phase jump a 
between them (here, we can fix 0 < a < n ), straightforward consideration of Eq. (1231) yields 

Weigen (L 0) = cos (a/2) — (T/2) sin a + O (i 2 ) , (27) 

the expansion of which at a —> 0 complies with Eq. USD. At a = 0, characteristic eigenfrequencies are 

Weigen (L = 0.5) » 0.90, W e igen (L = 1.0) « 0.74, W e ige„ (L = 1.5) M 0.61, W e igen (L = 2.0) w 0.51. (28) 

At large values of L , higher-order DMs appear, corresponding to n > 0 in Eq. (B^l) [higher-order linear DMs were 
also found numerically in the model including the defect profile © [44j], Indeed, higher-order DMs emerge, for each 
n > 1, at ui = 1, i.e., at 

=27171-a/2. (29) 

For instance, at a = 0 and with L slightly exceeding the first critical value E51) . i.e., L = 2-7T + l, 0 < l <C 27r, the 
emerging eigenfrequency is obtained as 4igen~ 1 — Z 2 /2, cf. Eq. (EH) , while its fundamental counterpart, found as 
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Fig. 1. Relations between the eigenfrequency of the fundamental defect mode, in the linear symmetric system with phase shift 
a [defined per Eq. lf3]l ]. and a, as produced by Eq. (1231) with n = 0, at fixed values of the width of the gapless layer: (a) 
L = 0.5; (b) L = 2.0. The dashed straight lines depict analytical approximation (1251) . 


a numerical solution of Eq. K51) . is w |;" gen °*(L = 2n) ps 0.22. As follows from Eq. (E51) . for given L the total number 
of defect modes is 


?*dm = 1 + 


L + (a/2) 
2n 


(30) 


where [...] denotes the integer part. 

Below, reporting numerical results for the full nonlinear system based on Eqs. ID and m , we focus on moderately 
wide gapless layers, with L < 2, and 0 < a < 7t/ 3, for which there is the single DM, corresponding to n = 0. 

Outside of bandgap Cl, an exact composite solution of the linearized equation (1^01) for delocalized states can be 
found too [cf. Eq. (12TT) ] : 


U( x) = e 


A + exp (iy/u> 2 — lx) + A_ exp (— i\fu> 2 — lx) , at x > 0, 
exp (iw ( x + -j)), at — L < x < 0, 

C + exp (iVw 2 — 1 (x + L )) + C- exp (— iy/u 2 — 1 ( x + L)) , at x < —L, 


where the amplitudes are 
A, 


C + = 


ei uL - e~i uL (w - Voj 2 - l) A 
2 Voj 2 — 1 (w + Voj 2 — l) 
e~i uL - e ^ L+ia (w - Vw 2 - l) 
2\/w 2 — 1 (w + Vco 2 — l) 


- \/w 2 - 1 j A* + , 

= (u> - \/a; 2 - l) e ia C* + . 


(31) 


(32) 


The existence of this solution at all w 2 > 1 corroborates that the insertion of the symmetric gapless layer does not 
make bandgap Cl broader than in the uniform grating, including the case when the phase shift is present, a/0. 

2.B.2. The asymmetric system 

The exact three-layer solution (ED-® for the linear DM can be easily generalized for the system with asymmetric 
reflectivity profile d): 


exp 


,«i 


W e 


- '/i-Weigen^ 


U(x) = eUi(x) = e ■ 


Bdm exp (iw e 


exp 


at x > 0, 
(x+ f)), at — L < x < 0, 


_1 (Tffir) + \! COs2 a - (Weigen ) 2 {x + L) 


(33) 


, at x < —L, 


S DM = exp 


2 '- C ° S 


(^eigen) 1 


■^eigen 


L ) 


(34) 
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and Eq. (l23l) for the eigenfrequency is replaced by 


^eie 


,L = 


1 T — 1 / \ . —if w eigen \ 

- COS (Weigen) + COS -— 

2 L V cos a J J 


2nn. 


(35) 


In the symmetric system, the DM exists at all values of width L of the gapless layer, including L = 0 [see Eq. 
urn provided that the phase shift satisfies condition a > —2 L, see Eq. (BUI) . The situation is drastically different in 
the asymmetric system, because, while the gapless layer attracts the DM (or soliton, see below), the reflectivity step, 
accounted for by cos a < 1 in Eq. repels the localized mode. The result of the competition of the attraction and 
repulsion is that, for given a, a solution of Eq. (1351) exists only if L exceeds a minimum value, 


T ^ T (DM) _ a + 47m 


(36) 


while at L < L m j n the DM does not exist. Exactly at L = L min , the eigenfrequency coincides with the edge of the 
respective bandgap m, Weigen = cosa. Further, it follows from Eq. EUH) that the total number of DMs in the 
asymmetric system is 


«DM 


1 + 


2Lcosa — a 


47T 


cf. the similar result (1301) for the symmetric model. Another corollary of Eq. (1551) is that, for fixed L , the fundamental 
DM (the one with n = 0) in the asymmetric system exists in interval 

0 < a < CKrnax, (37) 


with Qf m ax defined by equation 


2L COS (a max ) — CVmax ■ 


(38) 


At small and large L, solutions to Eq. (l38ll are 


f 2L-4L 3 , at L—> 0, 

\ (tt/2) [1 - 1/ (2 L )\, at L —> oo. 


(39) 


Particular solutions are 

amax (L = 0.5) « 0.74, a max (L = 1.0) w 1.03, a max (L = 1.5) w 1.17, a max (L = 2.0) m 1.25. (40) 

Finally, for given L , the DM’s eigenfrequency, Weigen, decreases from its value in the symmetric system [the one 
determined by Eq. (l23l) with a = 0] to the above-mentioned value, cos (a max (L)), with the increase of the asymmetry 
parameter, a , from zero to a max (L). These dependences are displayed in Fig. [2j which can be compared to their 
counterparts in the symmetric system with phase shift a, shown in Fig. [I] 

As mentioned above, in interval m between the broad and reduced bandgaps the asymmetric system supports 
semi-delocalized linear modes. The explicit form of these composite solutions is 

{ exp (| cos -1 oj — y/1 — u> 2 x ) , at x > 0, 
exp [t Q cos” 1 uj + uix)\ , at — L < x < 0, 

C+ exp (iy/uj 2 — cos 2 a (x + L )) + C- exp (— iy/u 2 — cos 2 a (x + L)) , at x < —L, 

where the amplitudes are 


~ (w + y/uj 2 — cos 2 a ) e l<?i — (cosa) e ^ 

= - / o -1 


a j 


— Vw 2 -' 


y/w 2 — i 


-C 


cos a 


+ > 


1 


= - cos 1 uj — Loj. 
2 


cf. Eqs. (1311) . (1551 and (1351) . 

















Fig. 2. (Color online) The relation between the eigenvalue of the fundamental linear defect mode in the asymmetric system 
and asymmetry parameter a, as per Eq. (1351) with n = 0, at fixed values of width L of the gapless layer: L = 0.5, L = 1.0, 
L = 1.5, and L = 2.0. The eigenvalues are shown in their existence region, 0 < a < a maX (T), see the text. The dashed green 
curve shows the boundary of bandgap m, u> = cos a, in the same region. 



2.C. The nonlinear correction to the eigenfrequency of the defect mode 

If the weak nonlinearity is taken into account, assuming that e in Eqs. (ED and ED is small but finite, in the first 
approximation it gives a small correction to the eigenfrequency, 


k^eigen 


o (0) 

eigen 


£ 


2 o/ 2) . 

eigen 


(41) 


We perform further analysis of this situation for the symmetric system, with a = 0. Multiplying the first equation of 
system © by solution ED. found in the linear approximation, and performing integration over dx yields, with regard 
to relation 0, 


( 2 ) 3 \Ui{%)\ 4 dx 

’ eiSen 2/_+~| u&tfdx 



(42) 


Here and in the subsequent numerical analysis, we used the Newton’s iteration method to construct nonlinear DMs, 
using the linear modes, given by Eqs. ED - ED. as the initial guess. The solution domain, with periodic boundary 
conditions, is \x\ < 60 (in all the cases, the width of localized modes is much smaller than the size of the domain). 
The resulting w e igen(£) curves for the nonlinear mode is plotted in Fig. [3] for two fixed values of L. The figure shows 
close agreement between the numerically found eigenfrequency and the analytical prediction produced by Eqs. ED 
and flSl) (the agreement is equally good at other values of L). 

Typical examples of the nonlinear DM, numerically found for w > 0 and w < 0, are shown in Figs. ©a,b,c). As 
demonstrated below [see Fig. ©a)], the nonlinear states, i.e., as a matter of fact, pinned gap solitons, displayed here 
in panels (a) and (b,c), are stable and unstable, respectively. In addition, Fig. 0d) displays an example of the linear 
DM. 


2.D. Gap solitons 

2.D.I. The symmetric system 

Stationary solutions of the nonlinear symmetric system ©, including phase shift a , for gap solitons pinned to the 
gapless layer, which are subject to restriction 0, can also be constructed as composite modes, by juxtaposing solutions 
(fill) in the semi-infinite gratings, and a plane-wave state in the gapless layer [cf. solution (12T1) for the DM]: 


U(x) 


3 cosh[(sin 6){x-\-^)—i0/ 2] 5 

-Bsoi exp 

/? e ia/2. 


i ^ |-B S oi| ^ (x + L) 


sin 6 


cosh[(sin 6)(x-\-L— £)— i6/2] ’ 


(43) 
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(a) (b) 

Fig. 3. (Color online) Chains of circles display numerically found eigenfrequencies of the nonlinear defect mode versus its 
amplitude in the symmetric system with a = 0, at different values of the width of the gapless layer: (a) L = 0.5; (b) L = 2.0. 
Blue solid lines represent the respective analytical approximation produced by Eqs. m and (S3). 




(a) (b) 




(c) (d) 


Fig. 4. (Color online) Typical examples of the nonlinear and linear defect modes in the symmetric system with a = 0 and 
L = 1.5: (a) w = 0.4016, e = 0.3; (b) w = -0.8715, e = 0.8; (c) w = 0, e = 0.5132; (d) lu = 0.6099, e = 0. Values of to 
corresponding to (a) and (b) are virtually exactly predicted by Eq. (1421) . The linear mode in (d) is produced by Eq. (1211) . 
Because the solutions obey relation 0, the curves depicting imaginary parts of U(x) and V(x) are identical, while their real 
parts have opposite signs. 
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Fig. 5. (Color online) An example of the intensity profile of a numerically found stable “flat-top” soliton, at a = 0, L = 1.125, 
uj = 0 . 1 . 


B soi — 


/£ to/2 sin 0 

V 3 cosh (£ sin 9 + *0/2) 


(44) 


where 9 is related to the soliton’s frequency by Eq. m, and £ > 0 is an offset, which is determined by the condition 
of the continuity of expression (l43l) at x = 0 [the continuity at x = —L is provided by Eq. (l44l) ]: 


cos 9 - 


2 sin 9 


cosh (2£ sin 9) cos 6 


L, + — = 






(45) 


with arbitrary integer n. The linearization limit corresponds to £ —> oo, which brings Eq. (1451) back into the form 
of Eq. (l23l) for the linear DM, taking relation (fl2l) into regard. This fact implies that the nonlinear extension of the 
linear DM coincides with the pinned gap solitons. In other words, it can be readily checked that, for given frequency 
uo belonging to bandgap m, the nonlinear DM coincides with the respective gap soliton. 

It is relevant to mention that exact composite solutions were previously obtained in other nonlinear models with 
multi-layer structures [55|-j57l|. Most of those models are built as nonlinear lattices of the Kronig-Penney type [58| . 

Thus, for given 9 = cos -1 (w so i), Eq. (Il5l) determines offset £ and, through it, the entire solution given by Eqs. 
and l| 44D for the fundamental pinned soliton at n = 0, and higher-order ones at n ^ 0. In particular, £ and 9 
determine the squared amplitude of the soliton: 


|£ sol | 2 


4 sin 2 9 

3 cosh (2£ sin 9) + cos 9 


(46) 


Although the analytical solution for the pinned soliton is found here in an implicit form, it produces some simple 
consequences. First, the intensity profile corresponding to Eq. (H51) is obviously a flat-top one, as \U{x)\ + |V(a;)| = 
2 |I3 so i| 2 does not depend on x in the gapless layer, see Eq. (147>1) and a typical example in Fig. [5] Second, it is obvious 
too that the solution obeys the symmetry condition m- 

Note that the composite solution given by Eq. (1431) is smoothest , i.e., without jumps of the derivative of the intensity 
profile, \U{x)\ 2 , at edges of the gapless layer, x = 0 and x = —L , if Eq. (l45l) yields £ = 0. It is easy to see that, for 
the fundamental pinned soliton (n = 0), this happens at 


uj = 2 + a/ (2L). (47) 

The soliton with this frequency may exist if w falls into bandgap m, i-e-, at a < —2L. Interestingly, this is precisely 
the case when the linear DM does not exist, pursuant to Eq. 

Lastly, in the limit of L = 0, i.e., for the junction between two semi-infinite gratings with phase shift a between 
them, it is easy to find an explicit solution of Eq. (USD: 


£ = (2sin0) 1 In 


/ tan {9/2) + tan (a/4) \ 
\tan {9/2) — tan (a/4 )) 


(48) 


which exists for 9 > a/2. Thus, solitons with internal frequency [see Eq. ([T21l ] w so i = cos 9 < cos (a/2) can be pinned 
to the interface between two BGs with phase jump a. The limit of w so i = cos (a/2) exactly corresponds to the linear 
mode pinned to the same interface, as per Eq. m- 
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(a) (b) 

Fig. 6. (Color online) Example of intensity profiles of numerically found solitons in the asymmetric system, (a) A stable soliton 
for a = 7t/4, L = 1.625, oj = 0.2. (b) An unstable soliton for a = 7t/3, L = 1.125, lo = —0.2 


2.D.2. The asymmetric system 

Composite flat-top solitons can be constructed for the asymmetric system too, following the pattern of ansatz (1431) . 
but it is then necessary to introduce two different offsets £ in the left and right semi-infinite gratings, which makes 
the composite analytical solution very cumbersome. However, the interaction of the BG soliton with the weak 
asymmetric interface, which is characterized by small a -C 7r/2 and L <tC sind in Eq. (U), can be analyzed by means of 
the simplest (adiabatic) version of the perturbation theory, which uses the Hamiltonian term (1101) , accounting for the 
Bragg reflectivity, to generate an effective potential of the interaction of a BG soliton with the interface, Wmt, while 
the perturbation of the soliton’s shape is disregarded [59]. A straightforward analysis yields the following asymptotic 
expression for the interaction force acting on the soliton with the center located at point 2 far to the left from the 
interface, —2sind 1: 


3 W- 8 

-Pint =-~ (sin 2 9) (4 L sin 9 - a 2 ) exp (-2 (sin 9) 2). (49) 

The soliton is captured by the weak interface if the the remote soliton is attracted to it, i.e., Fj nt > 0. Thus, it follows 
from Eq. flU) that the capture does not occur if the soliton’s amplitude is not large enough, sind < a 2 / (4 L). Because 
sin# cannot exceed 1, the final conclusion is that the interface cannot trap any BG soliton under condition 

L<L ( ^=a 2 / 4. (50) 

This prediction is qualitatively similar to that for the existence condition of the linear DM, given by Eq. m with 
n = 0, the difference being that, for a —> 0, ~ a 2 is much smaller than ~ a. 

3. Numerical results for gap solitons 

Solitons in the symmetric system (with the phase shift a) were obtained above in the implicit analytical form, based 
on Eqs. nm oni ' • Here, we focus on constructing a family of fundamental solitons [whose linear limit correspond to 
n = 0 in Eq. (1351) ] in the asymmetric system in a numerical form, by means of the Newton’s method. The results also 
include the symmetric system with a = 0 (the one with no phase shift). To this end, the exact gap soliton, given by 
Eq. (fill) for the uniform BG, was used as the initial guess. An essential conclusion of the numerical analysis, which 
complies with the similar property of the symmetric system [represented by the above-mentioned correspondence 
between Eqs. m and (E51) . as well as between (1151) and (EH)], is that nonlinear DMs, obtained as a continuation of 
their linear counterparts given by Eqs. El])-®- are tantamount to gap solitons pinned to the gapless layer. Stability 
of the so found solitons was then tested by means of direct simulations of their perturbed evolution, which were run 
up to t = 10000, using the split-step fast-Fourier-transform algorithm, with absorbers set at edges of the integration 
domain. 

Typical examples of numerically found stable and unstable composite solitons in the asymmetric system are pre¬ 
sented in Fig. [6] The asymmetry is clearly featured by shapes of the solitons, cf. the symmetric ones in Figs. [5] 
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(a) (b) 



CD 


(c) 


Fig. 7. (Color online) Energy E of the gap solitons versus u> in the asymmetric system for different L. (a) a = 0; (b) a = 7t/4; 
(c) a = 7t/3. Continuous and dotted segments denote stable and unstable solutions, respectively. In (a), the curve for L = 0 
represents the exact solution m, with E = (8/3) arccosw, as per Eqs. (1121) and (1141) . the boundary between the stable and 
unstable segments corresponding to Eq. m- In panels (b) and (c), the numerical solution cannot be extended beyond bottom 
points at which the branches terminate. 


Soliton families are characterized by the dependence of the total energy, E [see Eq. ©], on frequency w, which 
is shown in Fig. [7] for different fixed values of L. This figure also designates stable and unstable portions of the 
solution branches. Note that values of w at E — > 0 in the symmetric system with a = 0 [panel ©a)], which, obviously, 
correspond the linear DM, exactly agree with the respective eigenvalues (1751) . On the other hand, the fact that the 
solution branches cannot be extended up to E = 0, i.e., to the linear limit, at a > 0, and become progressively shorter 
with the increase of a and decrease of L, as seen in panels [71 b,c), is explained by the nonexistence of the linear DM 
in the asymmetric system for L too small, according to Eq. (l36l) . as well as by the nonexistence of the full nonlinear 
solutions at small L, according to Eq. (1501) . In terms of the numerical solution, the limitation on the existence of 
the localized modes imposed by the asymmetry is more severe than predicted by Eq. (1361) . because exactly at the 
existence limit given by this equation, L = a/ (2 cos a), the mode’s intrinsic frequency coincides with the edge of the 
bandgap, i.e., the mode is still delocalized and cannot be found in the numerical form in a finite domain. 

Figure [8] summarizes the numerical results by depicting stability boundaries for the solitons in the plane of ( L,E ) 
for different fixed values of a in Eq. ©. It was found that the existence and stability areas vanish at a > a CI ss n/2.8, 
which is explained by the above-mentioned trend to the disappearance of the localization modes with the increase of 
asymmetry a [obviously, no localized modes may exist at a = 7r/2, when bandgap (1131) with re = cos a shrinks to 
zero]. Further, in accordance with the above analytical predictions [see Eqs. m and ®], gap solitons cannot be 
pinned by the asymmetric layer whose thickness, L , is too small, while there is no such limitation for the symmetric 
layer, with a = 0. Finally, following the pattern of the gap solitons in the uniform BG [see Eqs. m and Cl], the 
pinned solitons tend to destabilize with the increase of the energy, and they do not exist at very large values of the 
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(a) (b) 



(c) 


Fig. 8. (Color online) Existence and stability regions for the pinned gap solitons in the plane of (L, E) of the asymmetric 
system at different fixed values of a: (a) a = 0 (the symmetric system); (b) a = 7r/4; (c) a = 7t/3. In the case of a = 0, the 
results for L = 0 correspond to Eqs. © and 09- 


energy. 

Evolution of unstable solitons was investigated by means of direct simulations. As shown in Fig. [9l unstable solitons 
are not destroyed. Instead, they shed off a part of their energy (which is eventually eliminated by the edge absorbers) 
and thus transform into stable solitons with smaller energy. 


4. Collisions of moving gap solitons with the gapless layer 

Collisions of moving BG solitons with the grating-free layer is another natural problem in the present setting. It was 
addressed by means of direct simulations, starting with a soliton initially set far from the layer, i.e., with the center 
placed at some large Xq > 0, and with initial velocity c < 0, so that it travels to the left, hitting the central layer. 
If the layer attracts the soliton, the collision may result in its splitting into reflected and transmitted pulses, plus a 
trapped mode. In the uniform grating with k = 1, the initial conditions for the moving gap soliton are taken as the 
exact solution (written here at the initial moment, t = 0) [&, 9j: 


^sol — 


^sol - 


2(1 + c) 
3- c 2 
2(1 —c) 
3 — c 2 


(l — c 2 ) W (X)exp [i(j>(X) — VT cos 9 ], 


(l — c 2 ) 1 ^ 4 W* (A)exp [i(j>(X) — iT cos 9] 


( 51 ) 
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Fig. 9. (Color online) A typical example of the evolution of an unstable soliton in the asymmetric system, at uj = —0.3, 
L = 1.125, a = 7t/4. The top panel shows the drop of the soliton’s total energy, due to radiation loss in the course of the 
transition to the stable shape. The eventually established soliton can be identified as one with u> = 0.2, which lies at the edge 
of the stability interval in Fig. Etc). 


X = (l — c 2 ) 1/2 x, T= -c(l-c 2 ) 1/2 t, 

cj>(X) = — —^-^tan -1 {tanh [(sin 9)X] tan (9/2)} , (52) 

o c 

W(X) = (sin 6) sech [(sin 9)X — i (9/2)] , 

the above solution m corresponding to c = 0. In this section we report results for the soliton with 9 = 0.29-7T, which 
is stable in the uniform BG. 

We first consider the layer with L = 1.0 and a = 0 (the symmetric system). Results of the collisions with this layer 
are shown in Fig. [10] A slow soliton, originally moving with velocity c = —0.1, splits into two parts, a reflected soliton 
and a trapped mode. The latter one keeps 23% of the total energy of the incident soliton. The incoming soliton with 
c = —0.5 splits into three parts: a reflected one, a trapped mode, and a weak transmitted radiation jet. In this case, 
the trapped mode keeps 13% of the total energy. Lastly, a fast incident soliton with c = —0.9 splits into two parts, a 
transmitted one, and a weak reflected radiation jet. 

As shown in Fig. El the difference in the case of the collision of the moving soliton with the gapless layer in the 
asymmetric system, with sufficiently large values of a, is that the trapped mode does not emerge as a result. This 
observation is explained by the reduction and disappearance of the stability region for the pinned solitons with the 
growth of a, as seen in Fig. [5] In particular, Fig. fill demonstrates that, at a = 7r/3, the slow incident soliton with 
c = —0.1 is reflected, and faster ones, with c = —0.5 and c = —0.9, are transmitted. In addition to these outcomes, 
radiation jets are generated too. Note that the transmission is facilitated, in this case, by the fact the incident soliton 
is moving from the domain with the stronger Bragg reflectivity to the one with the weaker reflectivity, according to 
Eq. ©. 

5. Conclusion 

We have introduced two models of the Kerr-nonlinear optical medium, which include two semi-infinite BGs (Bragg 
gratings), separated by the grating-free (gapless) layer of width L. In the symmetric model, the bandgap in both 
gratings is the same, while a phase shift between them, a, is possible. The other model is asymmetric, with different 
widths of the bandgaps on both sides of the gapless layer, due to different Bragg reflectivities . The model can be 
implemented in the temporal and spatial domains alike. The attraction of the fields to the gapless layer gives rise 
to a DM (defect mode) in the linear version of the systems, and pinned composite gap solitons in the nonlinear 
one. However, in the asymmetric system the attraction competes with repulsion from the reflectivity jump, therefore 
the DM and pinned solitons do not exist if L is too small in this system. Both the linear and nonlinear modes are 
constructed analytically and numerically in both systems, in exact or approximate forms. The stability region for 
the pinned solitons in the asymmetric system is identified by means of systematic simulations, demonstrating the 
shrinkage and disappearance with the increase of a, and disappearance at small L too. Collisions of freely moving 
solitons with the gapless layer were also studied, suggesting a possibility of the creation of the trapped mode, i.e., a 
pulse of standing light, as a result of the collision. 
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Fig. 10. (Color online) Collisions of moving gap solitons with the gapless layer in the symmetric system, with a = 0. (a) A 
slow gap soliton with velocity c = —0.1 splits into a reflected pulse and a trapped mode, (b) A gap soliton with c = —0.5 splits 
into a reflected pulse, a trapped mode, and radiation, (c) A fast gap soliton with c = —0.9 splits into a transmitted soliton and 
a radiation jet. 


The work can be naturally extended in other directions. In particular, while the present analysis was focused on 
the fundamental trapped modes, it may be interesting to consider higher-order ones, corresponding to n > 1 in Eqs. 
m, HMD, 63, ©, and 63, which exist if L is large enough. Another relevant extension is to address a similar 
model with other types of nonlinear terms, such as quadratic or resonant (two-level) ones. 
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